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Abstract. We briefly describe the current version of the PHOENIX code. We 
then present some resuhs on the modehng of Type II supernovae and show that 
fits to observations can be obtained, when account is taken for sphericaUy sym- 
metric, hne-blanketed, expanding atmospheres. We describe the SEAM method 
of obtaining distances to supernovae and briefly discuss its future prospects. 



1. The PHOENIX Code 



PHOENIX fsee iHauschildt &: Baroiil Il999l . and references therein) is a general- 
ized, stellar model atmosphere code for treating both static and moving at- 
mospheres. The goal of PHOENIX is to be both as general as possible so that 
essentially all astrophysical objects can be modeled with a single code, and to 
make as few approximations as possible. Approximations are inevitable (par- 
ticularly in atomic data where laboratory values for most quantities are un- 
known); however, the agreement of synthetic spectra with observations across a 
broad class of astroph ysical objects is a ve r y good consistency check. We have 
mode l ed Pl anets /BDs |Barman et al." '2002"; 'Allard et al."200r; Schweitzer et alj 
l2nnil . I2nn2l 'l. Cool stars (.Hauschildt et al...l999.: .Hauschildt et ah .1999 ). Hot 
Star s (/?CMa, eCMa, DeneblAufdenberg et al.lll998l . Il999l . I2nn2t ). «-Lvrae. No- 



Il997t). and all types of supernovae 



Lentz et al.l 120011 : iBaron et al.lll999l : 



vae ( Hauschildt et al.lll997l: [S chwarz et al 
fSNe la. Ib/c HP, H b lBaron et al...l995b< : 
iMitchell et ani2002l h 

PHO ENIX solves the radiative transfer problem by using the short-characteristic 
method (lOlson et al.lll987l :l(5i son fc Kunasdll987l 'l to obtain the formal solution 
of the special relativistic, spherically symmetric radiative transfer equation (SS- 
RTE) along its characteristic rays. The scattering problem is sol ved via the 
use of a band-diagonal approximation to the discretized A-oper ator ( Hauschildt! 
[1992; Olson & Kunasz 1987; Hauschil dt" Storzer. fc Baronlll994^ as our choice of 
the approximate A-operator. This method can be implemented very efficiently 
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to obtain an accurate solution of the SSRTE for continuum and line transfer 
problems using only modest amounts of computer resources. 

We emphasize that PHOENIX solves the radiative transfer problem including 
a full treatment of special relativistic radiative transfer in spherical geometry 
for all lines and continua. In addition we enforce the generalized condition of 
radiative equilibrium in the Lagrangian frame, including all velocity terms and 
deposition of energy from radiative decay or from external irradiation. 

We also include a full n on-LTE treatment of most ions, using model atoms 
constructed from the data o fE uruc 3 (|l99.\.1994a. .b) and/or from the CHIANTI 
and APED databases. The code uses Fortran-95 data structures to access the 
different databases and model atoms from either or both databases can be se- 
lected at execution time. 

Absorption and emission is treated assuming complete redistribution and 
detailed depth-dependent profiles for the lines. Fluorescence effects are included 
in the NLTE treatment. The equation of state used includes up to 26 ionization 
stages of 40 elements as well as up to 600 molecules. 



2. Early Supernova Spectra and Distances 

Early spectra of supernovae are important to model with synthetic spectra. 
Due to the geometrical dilution, early spectra tell us about the composition and 
velocity structure of the outermost layers of the supernova ejecta. Therefore, 
they can provide clues to the progenitor structure of SNe la, the primordial 
composition of all types of supernovae and the nearby circumstellar medium. 
iBaron et al. |(|2003) showed that detailed models of the early spectra of SN 1993J 
could reveal the primordial metalicity as well as the total reddening to the 
supernova. 

With all large computer codes, code improvement is an ongoing process. 
We recently improved the NLTE equation of state and Figure ^ compares the 
results with and without the improvement. The change is not dramatic and in 
the right direction, the new er model fits Ha better. 

We ( Baron et al.ll20do[ ) were also able to determine the reddening to SN 1999em, 



by modeling the early spectra. We found that E{B — V) = 0.1 it 0.05 for the 
total (foreground plus parent) color excess. 

In order to determine the values of the fundamental cosmological parame- 
ters, an accurate distance indicator visible to high redshift is required. Super- 
novae are extremely bright and hence can be detected at cosmological distances 
with modern large telescopes. The reliability of SNe la as distance indicators 
improved significantly with the realizatio n that the lum inosity at peak was cor- 
related with the width of the light curve (lPhilliDjll99l 'l and hence that SNe la 
were correctable candles in much the saine way that Cepheids are ()PhilliDs et al.l 
119991 : iGoldhaber et al.ll200ll : iBiess eranil99,'j^. This work and the development 
of highly efficient search strategies (|Perlmutter et al.lll997h sparked two groups 
to use SNe la to nieasure the deceleration param eter and to discover the dark 
energy ()Riess et al.lll998l : iPerlmutter et al.lll999l 'l. 

All of the work with SNe la is empirical, based on observed SNe la template 
light curves. Another method of determi ning distances using supernovae is the 
"expanding photosphere method" (EPM, iKirshner &: KwanI 19741 : iBranch et al.l 
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Figure 1. A comparison of synthetic spectra of SN 1993W on Aug 20, 1993 
using old (incorrect) and new (correct) NLTE equation of state. The change 
is not dramatic, but Ha is better fit with the new EOS. 
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Figure 2. The left panel displays a model with a lower temperature and 
extinction E{B — V) = 0.05, while the right panel shows a model with a higher 
temperature and extinction E{B — V) = 0.10. The Ca II H+K line is very 
temperature sensitive at this early phase and the higher extinction is clearly 
preferred. 



ll98ll : lEastman fc Kirshner"l989': ' Eastman et al.iri996r i a variation of the Baade- 
Wesselink method (■Baad&.1926. ). The EPM method assumes that for SNe IIP, 
with intact hydrogen envelopes, the spectrum is not far from that of a blackbody 
and hence the luminosity is approximately given hy L = An R"^ aT"^, where R 
is the radius of the photosphere, T is the effective temperature, a is the radiation 
constant, and C is the "dilution factor" which takes i nto account that in a scatter- 
ing do mi nated atmosphere the blackbody is diluted ( Hershkowitz. Linder. Sz Wagoned 
Il986al lbl: iHershkowitz fc Wagoned Il987j ) . The effective temperature is found 
from observed colors, so in fact is a color temperature and not an effective 
temperature, the photospheric velocity can be estimated from observed spectra 
using the velocities of the weakest lines, R = vt, the dilution factor is estimated 
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from synthetic spectral models, and t comes from the light curve and demanding 
self-consistency. 

This method suffers from uncertainties in determining the dilution factors, 
the difficulty of knowing which lines to use as velocity indicators, uncertainties 
between color temperatures and effective temperatures, and questions of how to 
match the photospheric radius used in the m odels to determine the dilution fac- 
tor a nd the radius of the line forming region ((Hamuv et al. l l2nnil : lLeon^ et al.l 
|2002). In spi te of this the EPM method was successfully applied to SN 1987A 
in the LMC ( Eastman Sz Kirshneilll989l : feranch 1987 ) which led to hopes that 
the EPM method would lead to accurate distances, independent of other astro- 
nomical calibrators. R ecently, the EPM m ethod wa s applied to the very well ob- 
served SN IIP 1999em ( Hamuv et al..i.200klLeonard et al...2002l:lElmhamdi et al.1 
I2OO3). Ah three groups found a distance of 7.5-8.0 Mdc. iLeonard et al.l (|2003l ) 
subsequently used HST to obtain a Cepheid distance to the parent galaxy of 
SN 1999em, NGC 1637, and found 11.7 ± 1.0 Mpc. 

With PHOENIX, accurate synthetic spectra of all types of supernovae can 
be calculated. The Spec t ral-fitting Expanding Atmosphere M ethod (SEAM) 
( Baron et al]ll995al . lT99^ : lLentz et ani200ll : iMitcheh et al.ll2002 l was developed 
using PHOENIX. While SEAM is similar to EPM in spirit, it avoids the use of 
dilution factors and color temperatures. Velocities are determined accurately by 
actually fitting synthetic and observed spectra. The radius is still determined 
by the relationship R = vt, (which is an excellent approximation because all 
supernovae quickly reach homologous expansion) and the explosion time is found 
by demanding self consistency. SEAM uses all the spectral information available 
in the observed spectra simultaneously which broadens the base of parameter 
determination. Since the spectral energy distribution is known completely from 
the calculated synthetic spectra, one may calculate the absolute magnitude, Mx, 
in any particular band X, Mx = — 2.5 log /q°° 5x(A) L>, dA + Cx, where Sx is 
the response of filter X, Lx is the luminosity per unit wavelength, and Cx is the 
zero point of filter X determined from standard stars. Then one immediately 
obtains a distance modulus ^ix = 'mx — Mx — Ax = 5 log (d/lOpc), where mx is 
the apparent magnitude in band X and Ax is the extinction due to dust along 
the line of sight either in the host galaxy or in our own galaxy. The SEAM 
method does not need to invoke a blackbody assumption or to calculate dilution 
factors. 



3. Fits and Conclusions 

Using the latest version of PHOENIX versi on 13.06, we have ca,lcula ted synthetic 
spectra of SN 1999em using Model S15 of lWooslev fc Weaver! (Il995l 'l. The model 
was expanded homologously and the gamma-ray deposition was paramet erized 
to be consistent with the nickel mixing found in SN 1987A (MitchelTeEaOEoOl). 
The abundances were taken directly from the model, and the effects of radioac- 
tive decay were taken into account. No attempt was made to adjust the met- 
alicity or helium mixing. The model had abundances: X = 0.64 Y = 0.34 
Z = 0.019. Figure ini shows the results are quite close to observed results except 
at the latest epoch where the blue is strongly dependent on the exact abun- 
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dances. With t hese models we are able to obtain a reasonably accurate distance 



Figure 3. The synthetic spectra are compared to observed spectra at 5 
different epoc hs. The observed spe ctra were obtained at CTIO for Oct 30, Nov 
2^and Nov 18 l|Hamuv et al.l200H) . at HST and FLWO on Nov 5 llBaron et al. 
[20 0011 and the optical spectrum on Nov 28 was obtained at Lick ({Leonard et al. 
I2OO2) while the IR was obtained at CTIO ijHamuv et alJl200H) . The observed 
fluxes have been offset for clarity. 



Results from both nearby and high-z searches require quantitative spec- 
troscopy to interpret their results. Results from hydrodynamical calculations 
also require quantitative spectroscopy to validate/falsify their results. The 
SEAM method applied to both nearby and high-z samples will yield interesting 
cosmological results. SEAM provides a route to obtaining cosmological param- 
eters completely independent from SNe la. 
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